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Using constraint-based local search, we effectively model and efficiently solve the problem of bal- 
ancing the traffic demands on portions of the European airspace while ensuring that their capacity 
constraints are satisfied. The traffic demand of a portion of airspace is the hourly number of flights 
planned to enter it, and its capacity is the upper bound on this number under which air-traffic con- 
trollers can work. Currently, the only form of demand-capacity balancing we allow is ground hold- 
ing, that is the changing of the take-off times of not yet airborne flights. Experiments with projected 
European flight plans of the year 2030 show that already this first form of demand-capacity balanc- 
ing is feasible without incurring too much total delay and that it can lead to a significantly better 
demand-capacity balance. 

1 Introduction 

The objective of air traffic management (ATM) is to ensure a safe, fair, and efficient flow of air traffic, 
under minimal environmental impact, subject to constraints on aircraft separation, airspace capacity, and 
airport capacity. The mission of EuroControl, the European Organisation for the Safety of Air Navi- 
gation (www.eurocontrol.int), is to promote the harmonisation of the different national ATM systems. 
EuroControl is the counterpart of the Federal Aviation Administration (R\A) of the USA. 

1.1 Air Traffic Management in Europe 

Current ATM systems in Europe are fragmented and already stretched to the limit, and hence unable 
to cope with the traffic volume foreseen for the year 2020 and beyond. Toward ensuring the required 
sustainability, an ambitious research effort at pan-European level and under the leadership of EuroControl 
is now underway to develop a fully integrated next-generation ATM system for the so-called Single 
European Sky ( www.sesarju.eu ). The context of this paper is our project on long-term innovative research 
with the EuroControl Experimental Centre (EEC) in Bretigny, France. 

Today, on the one hand, flight planning is made globally, at the strategic and tactical levels, for 
all 38 EuroControl countries by the Central Flow Management Unit (CFMU) of EuroControl, upon 
negotiation with the airlines but without sufficient effort at avoiding traffic bottleneck areas. On the 
other hand, flight control takes place more locally, within regional air-traffic control centres (ATCC), but 
without a super-regional view when flight re-planning has to be done. 

The operations of an ATCC rest upon a partition of its civilian airspace into sectors, that is three- 
dimensional (possibly concave) polygonal regions of airspace that are stacked at various altitudes (and 
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that often do not follow national boundaries). For each sector, a pair of air-traffic controllers (ATCo) 
try to ensure the sector capacity and aircraft separation constraints, at the operational level, and without 
necessarily physically being at the ATCC headquarters. For instance, the three BeNeLux countries and 
small bordering parts of Germany, France, and the North Sea are covered by the Maastricht ATCC in the 
Netherlands, whose airspace is partitioned into about 5 or 6 sectors at every layer, with 2 or 3 layers from 
the ground to upper airspace. 

Because of this fragmented mode of operation, the capacity of an ATCC is limited by its sector with 
the smallest capacity, the capacity of a sector being defined as the maximum hourly number of flights 
that may enter it. Indeed, the total capacity of an ATCC could be raised by an early identification of 
the traffic bottleneck areas and a redoing of the flight plans such that the traffic demand is more evenly 
balanced between its sectors. This has triggered the investigation of concepts dealing with multi-sector 
planning, where tools are developed to predict the traffic demand over several sectors and to manage the 
overall demand by anticipating peaks and proposing alternate plans. 

Furthermore, a long-term goal is to abandon the existing sectors and reorganise the entire airspace of 
the European Civil Aviation Conference (EC AC) by starting from a three-dimensional grid of same-sized 
box-shaped cells laid over the whole European continent and some nearby countries. For instance, in our 
experiments, using a grid of cells whose dimensions are 75 nm x 75 nm x 125 FLj^we have 4 layers 
with a total of 4600 ECAC cells, whose capacities are initially uniformly set to 40 entering flights per 
hour, to be enforced with a time step of 5 minutes. 

1.2 Dynamic Demand- Capacity Balancing 

This paper presents the first results on dynamic demand-capacity balancing at ECAC level, that is the 
dynamic modification of ECAC flight plans so as to satisfy the capacity constraints of all cells over a 
given time interval at minimal cost, thereby avoiding intolerable peaks of demand (and hence ATCo 
workload), as well as to balance the demands on all cells, thereby avoiding unacceptable dips of demand 
and unfair discrepancies between demands on the cells. There are many ways of modifying flight plans 
and defining the cost of such changes, and we postpone that discussion for a paragraph, so as not to 
disturb the presentation of the general objective. We are here only interested in aircraft that follow 
planned routes, rather than performing free flight. 

The tactical rolling-horizon scenario considered is as follows. At a given moment, suitably called 
now below, the (possibly human) demand-capacity manager queries the predicted traffic demands under 
the given flight plans for the cells of the entire ECAC airspace over a time interval that is some 30 
to 60 minutes long and starts some 3 hours after now. A look-ahead much below 3 hours would not 
give enough time for the implementation of some of the necessary re-planning. A look-ahead much 
beyond 3 hours would incur too much uncertainty in trajectory prediction (due to unpredictable weather 
conditions, say), and hence in demand prediction. If there are demand peaks, dips, or discrepancies that 
warrant interference, then the demand-capacity manager launches a re-planning process that suitably 
modifies the current flight plans. This process is to be repeated around the clock. For this to work, the 
time spent on re-planning should be very short, and the implementation effort of the modified plan should 
be offset by the resulting demand reductions and re-balancing among the cells. 

We now come to the announced discussion of how flight plans can be modified and how the cost of 
such changes can be defined. In this paper, we only consider ground holding, that is the changing of 
the take-off time of a not yet airborne flight, with the rest of its plan shifted forward accordingly. For 
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instance, the EEC suggested that our first experiments should allow ground holding by an integer amount 
of minutes within the range [0, . . . , 120]. Hence the total cost of such flight plan modifications could be 
the total delay incurred by all ground holding, and this cost is to be minimised. Our assumption here 
is that take-off times can be controlled with an accuracy of one minute, which is realistic nowadays. 
Under this choice, di flight plan for a given flight / is reduced to a time-ordered sequence of pairs c), 
indicating that / enters cell c at time Also, a modified flight plan is then obtained by adding a delay 
deltaT[f] to all the cell entry times of the original flight plan of /, hence there is a decision variable 
deltaT[f] G [0, . . . , 120] for every flight / of the original set of flight plans. Other flight plan modifications 
are discussed in the discussion of future work in Section IH 

In reality, there is another objective function, as the traffic demands on the cells should also be 
balanced over time and space. At the time of writing, it is not known yet how to combine any balancing 
cost with the ground holding cost into a single objective function, so we will just experiment with them 
separately. Quantifying the balancing cost may be done in many different ways, and we have decided to 
start with the standard deviation of the demands on the 4600 cells. 

In short, the optimisation problem of the present work can now be stated as follows: Given a set F of 
flight plans, a set C of airspace cells, and three moments now < ^ < ^ in time, return a modification (only 
by ground holding) of F of minimal cost, such that the traffic demands on the cells of C never exceed 
their capacities within the time interval [s^...^e\. In practice, an allocated amount maxlter of iterations 
is also given, and we want the best such flight plan modifications that can be computed within maxlter 
iterations. 



1.3 Contributions and Organisation of this Paper 

This paper makes the following main contributions. To our knowledge, this is the first time that multi- 
sector planning is attempted on the entire European airspace (and with even many more than the present 
20000 flights per day) and the first time that CBLS is applied to an ATM problem. We plan to make our 
constraints and objective functions progressively more realistic: see the discussion of related and future 
work in Section IH 

The rest of this paper is organised as follows. In Section [2} we summarise a constraint program 
modelling the dynamic demand-capacity balancing problem in a multi- sector-planning framework. In 
Section |3] we report on the experiments we made with that program. Finally, in Section |4j we conclude, 
discuss related work, and outline future work. 



2 Model 

Our constraint model is written in Comet |9|, an object-oriented constraint programming language with 
a back-end for constraint-based local search (available at www.dynadec.com). Some pre-processing 
(described below) via MySQL is made on the raw instance data, which consists of a possibly very large 
set of flight plans. 
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Figure 1: Sliding windows and the time-line 



2.1 Parameters and Notation 

The following notation is used to express the pre-processing and model, all times being expressed as 
integer amounts of minutes since some common origin: 



C = set of cells of some airspace 
F = set of flights 

df = planned departure / entry time of flight / in / into the airspace 
a/ = planned arrival / exit time of flight / in / from the airspace 
C 5 C/ = set of cells entered by flight / 
Cf^c = entry time of flight / into cell c 
now = time of launching the re-planning 
now < s = start time of interval for re-planning (typically 3 hours after now) 
s < e = end time of interval for re-planning (typically 30 to 60 minutes after s) 

w = length of time window over which metric is defined (60 minutes for capacity) 
/ = time step at which metric is checked (a divisor of ^ — ^, typically 5 minutes) 
m = {e — s)/t ^ number of time steps along the interval for re-planning 
g = maximum amount of ground holding per flight (currently 120 minutes) 
cap = capacity of a cell (40 entering flights per hour, currently for all cells) 



With the capacity constraint to be enforced every t minutes within the time interval . . . for re- 
planning, the number of flights entering a given cell has to be measured over a sliding window of w 
minutes that overlaps with [s^...^e\. Sliding window Sr (where < r < m) is the right-open time interval 
defined by (see Figure [T]): 

Sr — [s — w + r -t^. . . ^s + r -tl 
that is 5o = — w, . . . , ^[, 5i = — w + /,..., ^ + ^[, etc, until 5^ = — w, . . . , 
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2.2 Pre-Processing 

The pre-processing derives a few sets from the flight plans in order to express our constraint model and 
restrict search to only those flights and cells that are affected by the time interval [s^...^e] for re-planning: 

• Relevant Flights is the set of flights that are planned to depart / enter by the end e of the last sliding 
window and to arrive / exit upon the beginning ^ — w of the first sliding window, as otherwise they 
cannot be airborne in any sliding window of the interval for re-planning (see Figure [T] again): 

RF^if eF \ df<e & af>s-w} 

Note that the flights that can be ground-held so as to be airborne in some sliding window of the 
interval for re-planning need not be considered relevant, as they could only worsen the demands 
on the cells and as the objective of minimising the total amount of ground holding should set their 
ground holding to zero anyway. 

• Relevant Airborne Flights is the set of relevant flights that were planned to depart / enter by now 
and can thus not be ground-held: 

RAF ={feRF\df< now} 

• Relevant Waiting Flights is the set of relevant flights that are planned to depart / enter after now 
and are thus the only ones that can be ground-held: 

RWF ={feRF\df> now} 

• Sliding Window Scheduled Flights is the set of relevant waiting flights planned to contribute de- 
mand during a given sliding window Sr (where < r < m): 

SFr = {/ e RWF \ceCf & ef^c ^ Sr} 

• Sliding Window Potential Flights is the set of relevant waiting flights that can be ground-held in 
order to contribute demand during a given sliding window Sr (where < r < m): 

PFr^if eRWF\ceCf & s-w^r-t-g <ef^c<s-w^r-t} 

• Relevant Cells is the set of cells containing some possibly ground-held relevant waiting flight in 
some sliding window: 

RC={ceCf\feRWF & s-w-g<ef^c<e} 

Let F be a 2d array such that F[r, c] is initialised during the pre-processing phase to be the known part of 
the demand on relevant cell c during sliding window Sr, that is the number of relevant airborne flights 
that are inside c during Sr. The unknown part of that demand depends on the amount of ground holding 
imposed on the relevant waiting flights. 

2.3 The Decision Variables and the Initial Assignment 

Let delta! be a Id array of decision variables in the integer domain [0, . . . such that deltaT[f] is the 
time delay (duration of ground holding) for the relevant waiting flight /. All delays are set to in the 
initial assignment, as we want to ground-hold as few flights as possible. 
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Algorithm 1 Optimised constraint model 
f oralKr in 0. .m){ 
foralKc in RC){ 

currF = collect (f in cellFlight{c} : // currF = SF_r union PF_r 
f .enterCell >= (now + r*t) && f .enterCell < (s + r*t) 



if (P[r,c] + currF. getSizeO > cap) 
S.post (atmost (cap-P[r,c] , 
all(f in currF) 

(f .enterCell + delta! [f]) >= (s - w + r*t) && 
(f. enterCell + delta! [f]) < (s + r*t) 

)); 

} 

} 



2.4 The Constraints 

We currently only have the capacity constraints on all relevant cells during all sliding windows: 



For a sliding window Sr and a relevant cell c, this constraint takes a flight / from the union of the 
scheduled and potential relevant flights of Sr and reifies (for those who view this document in colour: by 
means of the red parentheses, which cast the truth of the embraced constraint into 1 and its falsity into 0) 
the interval membership to Sr of the actual entrance time (under ground holding) of / into c. The sum of 
the reified memberships for all such flights is the unknown demand on c during S/, augmented with the 
known demand P[r, c], the obtained total demand must never exceed the capacity (of the cell). 

Hence there would be as many as (m+ 1) • \RC\ constraints of the form ([T]), each involving \SFr U PFr\ 
decision variables and as many reified interval membership constraints for the sliding window Sr. How- 
ever, on the ECAC airspace, \RC\ easily goes into the thousands under the current values of the problem 
parameters. Fortunately, in practice, it is possible to remove about 85% of these constraints without 
losing any feasible solutions. Indeed, the constraints involving at most cap flights can be removed, since 
their total demand cannot exceed cap (recall that the total demand includes the unknown demand of the 
relevant waiting flights, which may or may not be ground-held sufficiently so as to enter the considered 
cell at the considered sliding window). In other words, the constraint ([T]) can safely be replaced by: 



) f; 



Vr (0 < r < m) . Vc G /^C. 

P[^^ c] ^HfeSFrUPFr i^f,c + deltaT[f] G Sr) < cap 



(1) 



thenP[ 



Vr (0 < r < m) . Vc G /^C. 
if P [r, c] + \SFrU PFr \ > cap 
y,c]+'£feSFrUPFr {ef^c + deltaT[f] E Sr) < cap 



(2) 



The model can further be shrunk by removing all the relevant cells that are only entered by relevant 
airborne flights. 

Our implementation uses the atmost constraint instead of reification, as seen in Algorithm [1} 
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2.5 The Objective Function 

The capacity constraints are considered soft during the search (in the sense that they may be violated by 
moves), hence the objective function must take both flight delays and constraint violations into account 
in order to ensure that the solution has minimum cost. The objective function to be minimised can be 
defined with weighted terms as follows: 

w • ^ deltaT[f] + v • violations {constmints ( 2 ) ) (3) 

feRWF 

where w and v are dynamically changing integer weights, both initially 1 and maintained by the search 
procedure (see Section |2.6| ) to minimise the sum of the delays. When the objective function cannot 
be improved, the values of the weights will be increased linearly. After finding the first solution, the 
diversification increases the violations and decreases the sum of the delays, so the balance between these 
weights must be reset to the initial values. 

In case we also want to re-balance the demands on the cells, we have to add a weighted term to the 
objective function ([3]), namely 

• {(^ceRC, o<r<mdemand[c,r]) 

where is an integer weight, demand[c, r] is the number of flights entering cell c during sliding window 
Sr under the currently chosen delays, and Oi^rKi) denotes the standard deviation of the expressions h{i) 
as / ranges over the set T. In first approximation, we will always have w • = 0, that is we will minimise 
only one of these terms of the objective function. 



2.6 The Search Procedure and Heuristics 

We here only give a heuristic and meta-heuristic for minimising the total ground holding (when w ^ 
& = 0). Different (meta-)heuristics may have to be designed for (simultaneous) demand balancing 
(when 7^ 0). 

The search procedure consists of a multi-state heuristic and a meta-heuristic. An event-based tabu 
search is used together with an exponential diversification in order to escape local minima during the 
constraint satisfaction phase, as well as an intensification in order to improve the objective function af- 
ter the satisfaction phase as the meta heuristic. We devised a heuristic with three states, and the search 
procedure may switch states while making moves. This heuristic is quite involved, but it gives bet- 
ter results than the other approaches we tried, including satisfaction with generic tabu search, variable 
neighbourhood search, variable-depth neighbourhood search, and weighted optimisation. Among these 
approaches, only satisfaction with generic tabu search was faster; however it could not minimise the 
objective function better than our three-state heuristic. Hence we now discuss only the latter. 

As mentioned in Section [2} an optimal solution must satisfy all the constraints and minimise the sum 
of all delays. Furthermore, an optimal solution is expected to have a maximum of flights with a zero 
delay, and the percentage of flights with a given delay is expected to decrease as the delay increases. 
Indeed, previous studies as well as discussions with CFMU and aircraft operators have shown that it is 
preferable to have many flights with short delays than a few flights with long delays of the same sum. In 
our heuristic, we make use of this feature to narrow down the neighbourhood by first selecting a delay 
d in the domain [0, . . . ,g] and then looking for the flight / that achieves the largest violation decrease if 
deltaT[f] is assigned d. The selection process of d cannot be random: it should follow the exponential 
probability to force the variables to fit best into the expected solution. 
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Algorithm 2 The multi-state heuristic 
float X = 1.3; 

float pr2[i in 1..12] = (x^^i) * (x-1 . 0) / (x^^ (13 . 0) -x) ; 
Closure step = closure { 
if (state == 1){ 

selectPr(i in 1 . . 12) (pr2 [i] ) 

select(d in ((12-i)*10) . . ((13-i)*10) : d>0) 

selectMin(f in RWF, ad = S . getAssignDelta(deltaT [f ] ,d) : 
tabu[f] <= it && S. violations (delta! [f ] ) > && 
delta! [f] != d && ad < 0)(ad){ 
delta! [f] := d; 

} 

} 

else if (state == 2){ 

selectMax(f in RWF: tabu[f] <= it && S. violations (delta! [f] ) > 0) 
(S . violations (delta! [f ] ) ) { 

selectMin(d in 0..g, ad = S .getAssignDelta(delta! [f ] ,d) : 
delta! [f] != d && ad < 0)(ad,d){ 
delta! [f] := d; 

} 

} 

} 

else if (state == 3){ 

selectMin(f in RWF, d in 0..g, 

ad = S.getAssignDelta(delta![f] ,d) : tabu[f] <= it && 
S. violations (delta! [f ] ) > && delta! [f] != d && ad < 0)(ad,d){ 
delta! [f] := d; 

} 

} 

>; 



We can define the required exponential probability distribution function by dividing the geometric 
series by the sum of its term: 

nx,y) = ""'J^'^} , with m<y<n& ^ f{x,y) = 1 (4) 

where m and n are the powers of the first and last terms in the geometric series, respectively; these values 
set the scale of the probability distribution. Variable x is a real number called the ratio of the series; 
it affects the speed of growth of the function. Variable 3; is an integer ranging from m to n. The delay 
domain [0, . . . ,g] being potentially large, we divide it into ten-minute-steps by making 3; range from 
m=ltOfz=12 when g = 120; the ratio was experimentally set to x = 1.3. Algorithm [2] depicts the 
Comet implementation of the heuristic. 

In the first state, when the violation of the constraints falls below an experimentally determined 
threshold of 300, the search procedure switches to the second state, which selects the most violated 
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variable and tries to minimise the violation by using the minimum value in its domain. This state reduces 
the violation of the constraints very quickly but disregards the total delay, whereas the first state aims 
at reducing the total delay. Finally, when the violation of the constraints falls below an experimentally 
determined threshold of 5, the third state is activated to expand the neighbourhood to all the variables, so 
that it is guaranteed that the next move is the best one toward the solution, as long as we are not at local 
minimum. The slow performance of the last state will not affect the heuristic since it only applies to the 
last few violated variables. 

The meta-heuristic employs diversification to escape local minima. The search procedure is depicted 
in Algorithm[3] The diversification is also moving toward a better objective function by selecting a set of 
flights and setting their delays to zero. The selection is not totally random as it uses a reverse sequence 
of probabilities derived from ^ with a ratio of 1.5. Once the constraints are satisfied, the diversification 
level is set to a higher value in order to ensure more modification to the flights' delays, and the state 
changes according to the new violations. This process continues to satisfy the constraints again, while 
the objective function is monitored to store the best solution so far before timing out or maximum number 
of iteration is reached. 

It is worth mentioning that, for a problem on such a large number of variables and constraints, 
achieving a smaller neighbourhood requires much more time per move, which is not acceptable in real- 
time conditions. Therefore, a simpler yet faster search procedure can reach better solutions. In our case, 
the algorithm can produce acceptable results in a very short time but it will not be able to converge to the 
optimal solution given enough time. 

3 Experiments 

The EEC has provided us with a generated data-set of foreseen traffic of the year 2030, namely approx- 
imately 50000 flights per day, including their entry times into the currently considered 4600 cells and 
1200 airports of the EC AC airspace. 

The experiments were performed with Comet (revision 2 beta) under Linux Ubuntu 9.0 (32 bit) on 
an Intel Core 2 Duo T7300 2.0GHz with 2MB cache and 4GB RAM (of which only 2GB are available 
to Comet). The code is allocating more than 1GB memory to solve the most complex instance. 

Our experimental results are based on two time intervals for re-planning. The first chosen interval is 
from ^ = 9 pm to ^ = 10 pm, which is a reasonably busy hour. The second chosen interval generates the 
most congested cells and corresponds to the flights departing from ^ = 5 pm to ^ = 6 pm; there are cells 
in this interval entered by more than 250 flights per hour. Table[T]shows statistics on the cell demands for 
these two intervals, taken every / = 12 minutes in a sliding window of w = 60 minutes. These statistics 
are measured before and after running the minimisation of the total delay, with now preceding ^ by 3 
hours, g = 120 minutes of maximum ground holding, = (no demand re-balancing in the objective 
function), and maxlter = 40000 iterations allocated; the minimum and median values are the same before 
and after optimisation. Note that all mean values have shrunk, except for the last sliding window of the 
first interval, because many flights were delayed to fall into that window. The last two columns confirm 
that the demands on the cells can indeed hugely exceed their capacity (maximum cap = 40 flight entries 
per hour), but also show that this capacity can actually be enforced by ground-holding flights for up to 
g = 120 minutes. 

The distribution of the delays among the relevant flights for the first interval is presented in Figure [2] 
as a histogram. Note the logarithmic scale of the vertical axis: the vast majority of the flights are not 
ground-held at all and, as wanted, the numbers of flights are more or less decreasing when the delays are 
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Algorithm 3 The search procedure and meta-heuristic 
while (it<maxlter) { 
call (step) ; 

if (oldViol == S. violations 0) steady++; 

else steady = 0; 

if (S. violations == 0){ 

maxDiverse=largeSteps ; 

state := 1; 

foralKf in RWF) tabu[f] = 0; 
if (obj < bestSolution){ 

bestSolution = obj ; 

sol = new Solution (m) ; 

} 

} 

else { 

maxDiverse = smallSteps; 

if (S. violations <= 5) state := 3; 

else if (S .violations <= 300) state := 2; 

} 

if (steady == diversif yLevel){ 
with atomic (m){ 

foralKi in . .maxDiverse) { 
selectPr(i in 1 . . 12) (pr2 [i] ){ 

select (f in RWF: delta! [f] > ((i-l)*10) && delta! [f] <= (i*10)){ 
delta! [f] := 0; 

»» 

steady = 0; 

} 

oldViol = S .violations ; 

} 

sol . restore ; 



increasing. 

Finally, Table |2] gives for both intervals the run-time, the actual number of iterations, the numbers 
of relevant waiting flights (that is the number of decision variables) and relevant airborne flights (which 
cannot be ground-held) at moment now, the total and average amounts of ground holding for all relevant 
flights, the change in standard deviation of the demands over all cells, and the initial number of constraint 
violations. It is very important to observe that, although the aim was only to minimise the total delay, 
there also is a 35% (respectively 43%) improvement in the standard deviation of the cell demands. About 
70% (respectively 50%) of the relevant waiting flights need not be ground-held at all; note that this 
percentage does not include the relevant airborne flights, nor the many flights that could possibly be 
ground-held so as to be airborne within the interval for re-planning. 

Among the parameters of our model (see Section [ZT] ), the maximum amount of ground holding per 
flight (currently set to g = 120 minutes) is the most interesting one to play with. Our first experiments 
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Table 1: Statistics on the ECAC cell demands before and after minimisation of the total delay on two 
chosen time intervals = 9 pm to ^ = 10 pm, and s ^ 5 pm to ^ = 6 pm) that are 3 hours ahead, 
with capacity enforced within . . . every ^ = 12 minutes over a window of w = 60 minutes, with a 
maximum g = 120 minutes of ground holding, and with = (no demand balancing in the objective 
function) 



Interval 


Run-time 


Iterations 


\RWF\ 


\RAF\ 


Total Delay 


Avg Delay 


Demand Dev 


Violations 


9-10 pm 


140 sec 


8,000 


4,295 


463 


65,457 min 


13.76 min 


-35% 


9,013 


5-06 pm 


620 sec 


22,000 


8,806 


811 


246,267 min 


25.61 min 


-43% 


45,985 



Table 2: Statistics for the two chosen time intervals: run-times, numbers of iterations, numbers of relevant 
waiting and airborne flights, total and average delays of all relevant flights, changes in the standard 
deviations of the cell demands, and initial numbers of constraint violations 

show that it can only be lowered to about g = 100 minutes without losing the satisfiability of the capacity 
constraints, though without significant deterioration of the run-time or total delay. 

Note that the trivial way to meet all the capacity constraints (namely by delaying a maximum of 
flights out of the time interval for re-planning, and thus exacerbating the situation for the next such 
interval) would be sub-optimal since we aim at minimising the total delay. 

4 Conclusion 

Summary. Constraint-based local search (CBLS) offers a very effective medium for modelling and 
efficiently solving the problem of optimally balancing the traffic demands on an airspace of adjacent 
cells, while ensuring that their capacity constraints are never violated. The demand on a cell is here 
defined as the hourly number of flights entering it. The currently only allowed form of re-planning is the 
changing of the take-off times of not yet airborne flights. Experiments with projected European flight 
plans of the year 2030 show that this can already lead to significant demand reductions and re-balancing 
at little cost. 

Ultimately, we hope to show that EuroControl can gain from a tight integration of planning and 
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Figure 2: Numbers of delayed relevant flights for each delay up to g = 120 minutes, bundled into 5- 
minute-intervals 



control. As discussed in the introduction, currently flight planning is made globally by its CFMU unit 
a long time in advance but without achieving optimal flow and under data estimates that are almost 
certainly incorrect, whereas flight control is made locally by its ATCC units when more precise data is 
available but without a global view, so that naive re-planning takes place. 

A manifestation of this lack of integration of planning and control is that the delays reported in 
Table |2] and Figure |2] would seem quite unacceptable (for instance, about 20% of the flights incur at 
least 30 minutes of delay). This is due to the poor cell demand distribution of the input flight plans (see 
Table [T]) and is thus an argument for the proposed tighter integration of planning and control. 

Related Work. Our previous project |[6l with EuroControl headquarters in Brussels (Belgium) is here 
extended from five airspace sectors (covering a bit more than the three BeNeLux countries) to the entire 
ECAC airspace, divided into box-shaped cells rather than the current sectorisation, and (initially) using 
capacity as a metric rather than some form of air-traffic complexity. Since this is a very large jump in the 
size of the input data, and since the technique to be developed even has to cope in real time with traffic 
volumes of the year 2030 (about 50000 flights per day, as opposed to the 20000 flights nowadays), our 
first idea was to use CBLS instead of the traditional constraint programming (by heuristic-guided global 
tree search with propagation at every node explored) that we deployed in the previous project. Since the 
chosen CBLS modelling language, namely Comet, also has constraint programming (by global search) 
and integer programming back-ends, we plan to keep the model and evaluate the actual necessity of this 
choice. 

The works closest to ours were tested on actual flight data of the year 1995 for all of France ||4l and 
the year 1999 for one French ATCC [2]. In further contrast to our present work (which was tested on 
projected flight data of the year 2030 for the entire ECAC airspace, divided into uniform cells rather than 
sectors), traditional constraint programming (rather than CBLS) was used toward satisfying the capacity 
constraints of all sectors, and the objective of the latter work was to minimise the maximum of all the 
imposed ground holdings (rather than their sum). 

Another important related work concerns the airspace of the USA The main differences with our 
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work are as follows. They have (at least initially) static lists of alternative routes to pick from for each 
flight and do not consider changing the time plans, whereas we currently only dynamically modify time 
plans. Their sector workload constraints limit the average number of flights in a sector over a given time 
interval (like our demand) but also the number of PIPs. No multi-sector planning is performed to reduce 
and re-balance the workloads of contiguous sectors. However, a notion of airline equity is introduced 
toward a collaborative decision-making process between the FAA and the airlines. 

Finally, there is related work on minimising costs when holding flights (on the ground or in the air), 
if not re-routing them, in the face of dynamically changing weather conditions |3 1. The main differences 
with our work are as follows. Their objective is cost reduction for airlines and airports, whereas our work 
is airspace oriented. They consider ground holding and air holding, whereas we currently only consider 
ground holding. Their dynamic re-routing is on the projected two-dimensional plane, whereas we do not 
re-route yet. Their sector workload constraints limit the number of flights in a sector at any given time 
(this is called sector load below), but there is no multi-sector planning. 

A lot of related work is about dealing with potentially interacting pairs (PIPs) of flights (see JB Ul 
for instance). However, this is an operational issue, whereas our work is at the tactical level, so we do 
not have to worry about the number of PIPs in the modified flight plans and hence we do not enforce any 
aircraft separation constraints. 

As this brief overview of related work shows, the whole problem of optimal airspace and airport 
usage by the airlines is very complex, and only facets thereof are being explored in each project. Our 
work intends to reveal some new facets, such as multi-sector planning. 

Future work. At present, this work is only a feasibility study, and this induced some simplifications 
in our experiments. For instance, the cells in our grid should have different capacities (rather than the 
uniform 40 we experimented with), just like the current sectors, which have roughly the same traffic 
volumes. Indeed, demand balancing under equal capacities does not make much sense, witness for 
instance low-level cells in areas without airports. However, our model can now be used to fine-tune the 
individual cell capacities. 

We must add side constraints to the current trunk of only capacity constraints in order to make things 
more realistic. For instance, since the re-planning happens for a time interval in the future, constraints 
will be needed to make sure no unacceptable traffic demand is generated before that interval. Additional 
constraints are also needed to make sure that the flight plan modifications can be implemented sufficiently 
quickly and cheaply, and that doing so is still offset by the resulting demand reductions and re-balancing 
among cells. For instance, the number of flights affected by the changes may have to be kept under a 
given threshold. This is where airlines and pilots will probably have to be factored in. 

Our constraint model does not enforce any notion of first-planned-first- served, which is a goal for 
CFMU and the FAA, since the delays would then become very huge. We will investigate situations where 
such a fairness notion can be applied nevertheless. 

Capacity (the maximum number of hourly flights that may enter a portion of airspace) is not the 
only metric we can enforce (under a given time step). Another common ATM metric is load, namely 
the maximum number of flights that can be simultaneously present in the portion of airspace. In a 
previous study [5 1, it was concluded that load constraints lead to less total delay than capacity constraints. 
Ultimately, a workload metric should be used (as in our |6|), estimating the traffic complexity as it is 
perceived by air- traffic controllers. 

Many other forms of flight re-planning can be imagined beyond the ground holding we here exper- 
imented with, such as the vertical re-routing of flights along the planned 2D route, the horizontal 2D 
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re-routing along alternative routes (from a list of fixed or dynamically calculated routes), accelerations 
and decelerations, and air holding. The cost function has to be adapted accordingly. Such additional 
forms of re-planning should only be introduced if they are warranted by additional cost gains and by 
computational feasibility. 

There is a lot of other future work to do before an early prototype like ours can be deployed in a 
tactical context. Its main current objective is therefore research strategic, namely to provide a platform 
where new metrics can readily be experimented with, and where constraints and objective functions can 
readily be changed or added. This motivated the choice of constraint programming as implementation 
technology, since the maintenance of constraint programs is simplified compared to ad hoc programs. 
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